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Abstract 

A sequential fourth generation of quarks and leptons is allowed by precision electroweak con- 
straints if the mass splitting between the heavy quarks is between 50 and 80 GeV. Although heavy 
quarks can be easily detected at the LHC, it is very difficult to detect a sequential heavy charged 
lepton, L, due to large backgrounds. Should the L mass be above 250 GeV, it can not be pair- 
produced at a 500 GeV ILC We calculate the cross section for the one-loop process e + e~ — > Lr. 
Although the cross section is small, it may be detectable. We also consider contributions from the 
two Higgs doublet model and the Randall- Sundrum model, in which case the cross section can be 
substantially higher. 
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I. INTRODUCTION 



One of the simplest extensions of the Standard Model is the addition of a new, sequential 
generation of fermions. Interest in this fourth generation has waxed and waned over the 
years. After precision measurements of the Z width showed that there are precisely three 
weakly interacting neutrinos [l| , it became clear that the neutrino mass of a fourth generation 
would have to exceed 45 GeV and interest faded. 

During the 90 's there was intensive study of the phenomenology of additional quarks 
and leptons which were not sequential . Many grand unified theories have additional 
fermions, such as vectorlike isosinglet quarks and leptons, additional vectorlike states arise 
in gauge-mediated supersymmetry breaking models, and many additional models contain 
mirrorlike fermions. These models are still of interest, but they do not require sequential 
fermions (although they can accommodate them). 

Interest in a sequential fourth generation faded further with studies of precision elec- 
troweak constraints. The recent Particle Data Group analysis claimed that "An extra gen- 
eration of ordinary fermions is excluded at the 99.999% CL on the basis of the S parameter 
alone" [l| ■ However, this analysis assumes a mass-degenerate fourth generation. Since one 
of the most striking features of the mass spectrum of the first three generations is the wide 
range of masses, such an assumption may not be justifiable. 

Analyses of the effects of a non-degenerate sequential fourthgeneration originally focused 
on the case where the neutrino mass was of O(50) GeV j^, E3, Q] and used 2001 electroweak 
data. A comprehensive analysis of the current status of precision electroweak fits and a 
fourth generation was recently carried out by Kribs, Plehn, Spannowsky and Tait 8j. They 
noted that the constraints on the oblique parameters from combined electroweak data have 
been determined by both the LEP Electroweak Working Group js| and the Particle Data 
Group The two groups used somewhat different datasets and differ by roughly one 
standard deviation (see Ref. 8[ for a detailed discussion of the differences). Kribs et al. 
used the LEP Electroweak Working Group results, and found that a substantial region of 
fourth-generation parameter space is in agreement with all experimental constraints. In this 
region of parameter space, the mass splitting between the U and D quarks is between 50 
and 80 GeV. Bounds on the mass splitting between the charged lepton, L, and the neutrino, 
N, are less constrained since one considers both Dirac and Majorana neutrino masses [lol |. 
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Thus, we see that precision electroweak data do not exclude a sequential fourth gener- 
ation. U and D quark production at the Large Hadron Collider (LHC) will be relatively 
easy to detect. However, the heavy charged lepton, L, will be substantially more difficult to 
detect, primarily due to large backgrounds. Early LHC and SSC studies [ill . made the 
assumption of a massless fourth generation neutrino, and still concluded that detecting a 
heavy lepton with a mass greater than 250 GeV might not be possible. Calculations of heavy 



lepton production (with a heavy neutrino) exist [13|, but do not include any discussion of 
signatures or backgrounds. Therefore, it is likely that a heavy charged lepton with a mass 
greater than 250 GeV will not be detected at the LHC. 

At the International Linear Collider (ILC), heavy leptons can be easily produced and 
detected up to the kinematic limit. However, the initial stage of the ILC will probably 
be at a center of mass energy of 500 GeV, in which case pair production of heavy leptons 
with masses above 250 GeV will not be possible. The only possible production mechanism 
would be through single L production, in association with a lighter Standard Model charged 
lepton. Since mixing between L and /i or e is expected to be small, we will focus on the 
process e + e~ — > Lr, which can occur through a nonzero #34 mixing angle. Although single 
production of heavy charged leptons has been studied before , all of these studies 
considered vectorlike or mirrorlike leptons, and we know of no calculations of this process 
with a heavy neutrino at a linear collider. An analysis of sequential heavy charged leptons 
in Z-decays jlj] ignored the mass of the heavy neutrino. 

In the next section, we present the relevant diagrams in the sequential standard model 
and calculate the cross section as a function of the N and L masses. Section 3 considers the 
cross section in the two Higgs doublet model and in the Randall- Sundrum model. Finally, 
in Section 4, we discuss detection possibilities and present our conclusions. 

II. SEQUENTIAL STANDARD MODEL 

A single charged heavy lepton can only be produced if the fourth generation mixes with 
the lighter generations. Bounds on the mixing angle #34 arise from observation of universality 
in r decays; a nonzero mixing angle would multiple the rate by cos 2 #34. This was analyzed 



by Swain and Taylor 



151 ] who found a model-independent bound of sin 2 #34 < 0.007. This is 



a particularly interesting value. If one has a Fritzsch-type 2x2 mass matrix for the r and 
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L, one might expect sin 2 #34 to be approximately m r /m^, which gives 0.007 for an L mass 
of 250 GeV. We will assume this value of the mixing angle in our numerical results, and can 
easily scale the cross section for smaller mixing angles. 




SM. The Feynman gauge is assumed and the light electron-Higgs couplings are neglected. 

The diagrams are listed in Figures 1-3 and grouped as the self energy, vertex and box 
type contributions, respectively. We use the 't-Hooft-Feynman gauge throughout, and thus 
charged Goldstone bosons, G, must be included. Note that the electron-Higgs couplings are 
neglected due to small Yukawa couplings. The internal neutrino lines get a contribution 
from each of the four neutrinos, and thus each diagram is proportional to V^V^. When 
summing over the four neutrinos, parts of the matrix elements that are independent of 
the neutrino mass will cancel by unitarity of the 4-D CKM-like matrix. This causes the 
ultraviolet divergences to cancel in the sum over neutrinos. 

The calculation of the cross section is performed by using the FeynArts, FormCalc, and 
LoopTools packages [l?}]. We first patched the SM and 2HDM model files of the FeynArts 
package by introducing the fourth generation leptons and their interactions. Then, the 
numerical analysis was carried out in Fortran with the help of FormCalc and LoopTools. 
The cancellation of the ultraviolet and infrared divergences has been checked numerically 
and the expected cancellation was confirmed. In addition, as a separate check the expected 
null result for the cross section due to unitarity of the mixing matrix was also tested 
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FIG. 2: The leading order vertex contributions to the e + e — > Lt process in the sequential SM.The 
't-Hooft-Feynman gauge is assumed and the light electron-Higgs couplings are neglected. 




FIG. 3: The leading order box contributions to the e + e — > Lt process in the sequential SM. The 
't-Hooft-Feynman gauge is assumed and the light electron-Higgs couplings are neglected. 

numerically by setting the heavy neutrino mass to zero. Note that the same technique 
is applied for the calculation in the 2HDM, presented in the next section. 

The results are plotted in Figure 4 for neutrino masses of 300, 400 and 500 GeV. We see 
that, for sin 2 # 34 = 0.007, cross sections of the order of a few attobarns can be expected. We 
can also show that as the neutrino mass increases, the cross section grows rapidly, reaching 
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FIG. 4: The total cross section of e + e~ — > Lf as a function of the heavy lepton mass mj, for 
y/s = 500 GeV and various heavy neutrino masses in an unpolarized electron-positron beam within 
the Sequential SM framework. 

500 attobarns at = 2000 GeV. This is not surprising since the theory is chiral. Of 
course, the cross section scales with V^; the value we have chosen is the maximum allowed 
from the analysis of Swain and Taylor. 

The structure of the curves in Figure 4 can be easily understood. Since the theory is 
chiral, one expects the cross section to increase as the mass of the heavy neutrino 
increases. However, as seen from Figure 4 that this is not necessarily true for neutrinos in 
the 300 — 400 GeV mass range. One can understand, for example, why the the curve for 
ttitv = 300 GeV crosses and becomes bigger than the one for = 400 GeV and similar 
behavior occurs between the = 400 GeV and = 500 GeV curves. This is simply 
due to the fact that both the W boson and the heavy neutrino N go on-shell in the loop 
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if the condition rriL > "irtw + m N is kinematically satisfied. When is large enough to 
produce the W and N on-shell, the loop integrals develop imaginary parts, which can be 
calculated by using the Cutkosky rules, and results in enhancement of the cross section. 
One can calculate this by cutting through the W boson-z/j propagators (for i — 4) at the 
heavy lepton's leg in Figures 1-3. Thus, for example, the peak due this enhancement for the 
m>N = 300 GeV curve occurs at around + m\y and it shifts to the right for the tun = 400 
GeV curve. 

Are these small cross sections detectable? With an integrated luminosity of an inverse 
attobarn, expected at the ILC's full luminosity for a couple of years, one expects a handful of 
events. The tau is monochromatic, and is opposite a monochromatic W and a light neutrino. 
We know of no backgrounds to this signature, and a complete analysis would be worthwhile. 



III. THE TWO-HIGGS DOUBLET AND RANDALL- SUNDRUM MODELS 



A. The Two-Higgs Doublet Model 

The minimal Standard Model Higgs sector consists of one complex Higgs doublet. One of 
the simplest and most popular extensions of the Higgs sector is the two-Higgs doublet model 
(2HDM). By requiring that all fermions of a given electric charge couple to no more than 



one Higgs doublet [la], one avoids flavor changing neutral currents. This is accomplished 
with a simple Z 2 symmetry. 

The 2HDM is an attractive model for several reasons: 

• it contains charged Higgs bosons and pseudoscalars 

• it adds relatively few new arbitrary parameters 

• it allows for spontaneous CP violation, and can give sufficient baryogenesis 

• this structure of the Higgs sector is required in low-energy supersymmetric models 



A very detailed discussion of the 2HDM can be found in the Higgs Hunter's Guide 

This model has two complex, Y — 1, SU(2)i doublet scalar fields <f>x and <p2- The 
vacuum expectation values of the neutral components of the Higgs doublets are v\ and t>2, 
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respectively. It is useful to define 

tan/3 = — . (3.1) 

Vi 

The physical Higgs fields consist of two neutral scalars, a neutral pseudoscalar and a charged 
Higgs scalar. In the charged sector, there will be both a Goldstone boson and a physical 
Higgs state. The charged Higgs is given by 

E ± = -<^sin/3 + 0±sin/? (3.2) 




12 3 4 

FIG. 5: The extra self energy diagrams contributing to e + e~ — > Lf in 2HDM. 
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FIG. 6: The extra vertex diagrams contributing to e + e~ — > Lf in 2HDM. 



For our calculation, the neutral scalars will not contribute. However, the charged Higgs 
boson will contribute. One simply replaces the charged Goldstone boson G in Figures 1-3 
with the charged Higgs boson; these diagrams are shown in Figures 5 and 6. The only 
exception is that the ZW T H ± vertex vanishes [3]. There are now two new parameters in 
the calculation, the mass of the charged Higgs boson and tan (3. 

There are two versions of the 2HDM. In Model I, all of the fermions couple to one of 
the Higgs doublets; in Model II (which is included in supersymmetric models), the neutral 
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leptons couple to one doublet and the charged leptons couple to the other. The relevant 
Yukawa couplings are 



L-N -H + 
d-N -H + 
L-Ui- H + 



if 



2\/2mw sin 9w 

2y / 2m w sin W 
ie5 i)3 Vu 



tan /5 
m N 
tan /5 

m L y(i + 75 ), 



(l- 75 ) + m L F(l + 75 ) 
(I-75) +m ei Y(l + 7 5 ) 



(3.3) 

ly 2mw sin Uw 

where Y — — 1/ tan/5 for Model I and tan/5 for Model II and the vertices for the ordinary 
lepton - ordinary neutrino - H + can be found in the Higgs Hunter's Guide 181 ] . 

Constraints from b — > S7 force the mass of the charged Higgs to exceed approximately 



200 GeV 19|. tan/3 and cot/9 must be less than about 3 so that the charged and neutral 
lepton Yukawa coupling remain perturbative. 
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FIG. 7: The total cross section of e + e~ — > Lf as a function of the heavy lepton mass mi for 
^Js = 500 GeV in an unpolarized electron-positron beam for various tan p values in 2HDM. In both 
graphs, the heavy neutrino and charged Higgs masses are set 400 GeV and 200 GeV, respectively. 

The results are presented for Model I and for Model II in Figure [3 We see small changes 
in the cross section for 1 < tan/5 < 3, but substantial changes for 1 < cot/5 < 3. In both 
models, the cross section can be enhanced by up to a factor of ten, leading to much easier 
detection at the ILC Note that the vertices involving the heavy neutrino scales a 1/ tan/3, 
and thus the cross section is enhanced if tan /5 < 1 as seen from Figure [71 
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B. The Randall-Sundrum Model 



Another popular extension of the Standard Model is the Randall-Sundrum (RSI) model 

n 

[201 ] . In this model a slice of AdS^ space with curvature k is compactified on an Si/Z 2 
orbifold with radius R. The fixed points of the orbifold correspond to two 3-branes with 
opposite tensions. The 4D graviton is localized on the positive tension, or Planck brane, 
while the Higgs field is localized on the negative tension, or TeV brane. This model generates 
a gauge hierarchy given by the warp factor, e~ nkR , and for kR ~ 12 completely solves the 
hierarchy problem (see 21[ for a bound on kR from the Casimir force measurements). 

In the original formulation, the model had all of the fermions on the TeV brane. More 
interesting phenomenology can occur when the Standard Model fermions and gauge bosons 
can propagate in the bulk 22J. In this case, the profiles of bulk fermion wavefunctions 
depend on their 5D mass parameters. By choosing the lighter fermions to live near the 
Planck brane, one can naturally explain the small Yukawa couplings for the light fermions, 
since their overlap with the TeV-brane localized Higgs boson is exponentially supressed. 
Thus the model can also explain the flavor hierarchy, since large differences in Yukawa 
couplings can arise from small differences in the mass parameters. The flavor hierarchy 
simply becomes a matter of geography in the fifth dimension. 



In an interesting series of papers, Agashe, Perez and Soni [23|, [24|, [25] discussed the 
phenomenological implications of the flavor structure of these models. They noted that 
one expects larger flavor changing neutral currents for the heavier generations, thus evading 



bounds involving light quarks. In particular [25j, Agashe, et al. considered top flavor 
violation at colliders, considering t — > cZ at the LHC, and e + e~ — > tc at the ILC. Clearly, 
a similar process could lead to e + e~ — > Lt as well, if a fourth generation exists. The 
mechanism is caused by the fact that the couplings of the fermions to the gauge boson 
Kaluza-Klein (KK) modes are not universal due to the different profiles for the fermions, 
and mixing between the gauge KK modes and the gauge bosons leads to flavor violating 
couplings of the Z. We refer the reader to Ref. 25J for details. 



One can simply carry over the calculation of e + e~ — > tc in Ref. [25J to this model. There 
is, however, one crucial difference. Since the couplings of the left-handed b quark to the Z 
are measured to an accuracy of less than one percent and the b quark is in a doublet with the 
left-handed top, one can not put the left-handed top and bottom too close to the TeV brane. 
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The right-handed top, however, can be close to the TeV brane. Thus the top flavor violation 
is predominantly right-handed. In the four generation case, there are no such restrictions, 
therefore the L flavor violation is relatively unconstrained. For definiteness, we choose the 
same magnitude for the left- and right-handed flavor violation, and set the coefficient of the 
I/JutZ' 1 term to be the same as that of the t^^c^Z^ term in the Agashe, et al. analysis 



251 ] . This is not unreasonable, since m c /m t ~ m T jm^ indicates that similar mixing angles 



may be expected. 

Using this flavor violating coupling, one can find the total cross section for e + e~ — > 
Lt. The result depends as well on the KK scale. It has been shown 26| that a custodial 
SU(2) symmetry in the bulk can allow the KK gauge boson mass to be as low as 3 TeV, 
and perhaps somewhat lower if a modest fine-tuning is allowed, without conflicting with 
precision electroweak results. Rather than calculate the interference with the Standard 
Model diagrams, we simply look at the RS model effects in isolation. This is because the 
uncertainty in the flavor-violating couplings preclude precise calculations. We find that if 
Mkk is 1 TeV, then the cross section varies from 1.0 to 0.5 femtobarns as the L mass varies 
from 250 GeV to 350 GeV, and scales as 1/M^ K . Thus, we see a significant enhancement of 
the cross section in the KK mass range of 1 — 3 TeV. One should keep in mind that the KK 
gauge bosons, if they exist, will be discovered at the LHC long before the ILC is constructed. 



IV. DETECTION AND CONCLUSIONS 



There are two possible decay modes for the L. It can decay into NW, or into u T W . 
Of course, if the N is heavier only the latter decay is possible. Regardless, there will be 
substantial missing energy in the decay. 

For the e + e~ — > Lt process detection should be extremely straightforward, since the r 
is monochromatic. For an L mass of 300 GeV, the r energy is 160 GeV, leading to a decay 
length, 7ct, of 0.8 centimeters. This is comparable to the size of the inner vertex detector 
at the ILC. 

In a wide region of the mass-mixing angle plane, the L will decay into a v T and a W. 
This would seem to give a clear signature, with a monochromatic r, a monochromatic W 
and missing energy. The primary background will be from r pairs, where one of the r's is 
misidentified. A detailed Monte Carlo analysis is beyond the scope of this paper, but if the 
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background can be eliminated, then a few events will suffice to discover the L. 

In the Standard Model case, we have found that there will be a few events produced at 
the ILC, and the question of whether or not the L can be detected depends on the details 
of the detector and Monte Carlo simulations. 

We then considered contributions from the charged Higgs boson of a 2HDM, as well as 
flavor changing effects in the Randall-Sundrum model. In both cases, there are regions of 
parameter space in which the cross section is substantially higher, leading to straightforward 
detection at the ILC. 

Long before the ILC is built, the LHC will have determined whether or not a fourth 
generation exists. If it does exist, then detection of the charged heavy lepton at the LHC 
will be very difficult and perhaps impossible. At the ILC, if the mass of the heavy lepton 
is more than half ^/s, pair production will be impossible, and the process calculated in this 
paper may be the only mechanism for detection. 
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National Science Foundation, Grant PHY-0554854. The work of I.T. is supported by the 
NSERC of Canada under the Grant No. SAP01105354. 
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